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Two fundamental ingredients play a decisive role in the foundation of fluctuation relations: the 
principle of microreversibility and the fact that thermal equilibrium is described by the Gibbs 
canonical ensemble. Building on these two pillars we guide the reader through a self-contained 
exposition of the theory and applications of quantum fluctuation relations. These are exact results 
that constitute the fulcrum of the recent development of nonequilibrium thermodynamics beyond 
the linear response regime. The material is organized in a way that emphasizes the historical 
connection between quantum fluctuation relations and (non)-linear response theory. We also 
attempt to clarify a number of fundamental issues which were not completely settled in the prior 
literature. The main focus is on (i) work fluctuation relations for transiently driven closed or 
open quantum systems, and (ii) on fluctuation relations for heat and matter exchange in quantum 
transport settings. Recently performed and proposed experimental applications are presented and 
discussed. 
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I. INTRODUCTION 

This Colloquium focuses on fluctuation relations and 
in particular on their quantum versions. These relations 
constitute a research topic that recently has attracted a 
great deal of attention. At the microscopic level, matter 
is in a permanent state of agitation; consequently many 
physical quantities of interest continuously undergo ran- 
dom fluctuations. The purpose of statistical mechanics is 
the characterization of the statistical properties of those 
fluctuating quantities from the known laws of classical 
and quantum physics that govern the dynamics of the 
constituents of matter. A paradigmatic example is the 
Maxwell distribution of velocities in a rarefied gas at equi- 
librium, which follows from the sole assumptions that 
the micro-dynamics are Hamiltonian, and that the very 
many system constituents interact via negligible, short 
range forces (Khinchin 



1949). Besides the fluctuation 
of velocity (or energy) at equilibrium, one might be in- 
terested in the properties of other fluctuating quantities, 
e.g. heat and work, characterizing non-equilibrium trans- 
formations. Imposed by the reversibility of microscopic 
dynamical laws, the fluctuation relations put severe re- 
strictions on the form that the probability density func- 
tion (pdf) of the considered non-equilibrium fluctuating 
quantities may assume. Fluctuation relations are typi- 
cally expressed in the form 



Pf{x) = Pb{~x) exp[a(a; - b)], 



(1) 
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FIG. 1 Example of statistics obeying the fluctuation relation, 
Eq. |T|. Left panel: Probability distribution p(q) of number 
q of electrons, transported through a nano-j unction subject 
to an electrical potential difference. Right panel: the lin- 
ear behavior of \n[p{q) /p{—q)] evidences that p{q) obeys the 
fluctuation relation, Eq. (ij). In this example forward and 
backward protocols coincide yielding pb = Pf = P, and con- 
sequently 6 = in Eq. ([Tl. Data courtesy of Utsumi et al 
(|20T0f. 



where pf{x) is the probability density function (pdf) of 
the fluctuating quantity x during a nonequilibrium ther- 
modynamic transformation - referred to for simphcity 
as the forward {F) transformation -, and pb{x) is the 
pdf of X during the reversed (backward, B) transforma- 
tion. The precise meaning of these expressions will be 
amply clarified below. The real-valued constants a, &, 
contain information about the equilibrium starting points 
of the B and F transformations. Figure [l] depicts a prob- 
ability distribution satisfying the fluctuation relation, as 
measured in a recent experiment of electron transport 



through a nano-junction (Utsumi et al. 2010). We shall 



analyze this experiment in detail in Sec. VI 

As often happens in science, the historical development 
of theories is quite tortuous. Fluctuation relations are no 
exception in this respect. Without any intention of pro- 
viding a thorough and complete historical account, we 
will mention below a few milestones that, in our view, 
mark crucial steps in the historical development of quan- 
tum fluctuation relations. The beginning of the story 
might be traced back to the early years of the last cen- 



stein 


( 


1905 


1906a|bl 


Nyquist 


( 


1928 


) later, 



tury, with the work of Sutherland ( 1902 1905) and Ein- 



le linear 

response of a system in thermal equilibrium as it is driven 
out of equilibrium by an external force, is determined 
by the fluctuation properties of the system in the initial 
equi librium s t ate. Specificall y, Sutherland (1902 1905) 
and Einstein (1905 |1906a|b ) found a relation between 
the mobility of a Brownian particle (encoding informa- 
tion about its response to an externally applied force) 
and its diffusion constant ( encoding informa tion about its 
equilib rium fluctuations). Johnson] ( [192 8*) and Nyquist' 
(1928)^ discovered the corresponding relation between 



^ Nyquist already discusses in his Eq. (8) a precursor of the quan- 
tum fluctuation-dissipation theorem as developed later by lCallen] 
I and Welton| ( |1951^ . He only missed the correct form by omit- 
ting in his result the zero-point energy contribution, see also in 



the resistance of a circuit and the spontaneous current 
fluctuations occurring in absence of an applied electric 
potential. 



The next prominent step was taken by C alien and Wel- 



ton| ( 1951 ) who derived the previous results within a gen- 
eral quantum mechanical setting. The starting point of 
their analysis is a quantum mechanical system described 
by a Hamiltonian 'Ho- Initially this system stays in 
a thermal equilibrium state at the inverse temperature 
j3 = {kBT)~^ , wherein ks is the Boltzmann constant. 
This state is described by a density matrix of canonical 
form; i.e., it is given by a Gibbs state 



-mo 



(2) 



where Zq — Tr e~^'^° denotes the partition function of 
the unperturbed system and Tr denotes trace over its 
Hilbert space. At later times t > Q the system is per- 
turbed by the action of an external, in general time- 
dependent force At that couples to an observable Q of 
the system. The dynamics of the system then is gov- 
erned by the modified, time-dependent Hamiltonian 



■H{\t) = nQ-\tQ- 



(3) 



The approach of Callen and Welton (1951) was further 



systematized by Green ( 1952 1954) and in particular by 



Kubo (1957) who proved that the linear response is de- 
termined by a response function 4>BQ{t)i which gives the 
deviation {AB{t)) of the expectation value of an observ- 
able B to its unperturbed value as 



(ABit)) = [ (/.BQ(t -s)A,ds. 

J —oo 



(4) 



Kubo (1957) showed that the response function can be 
expressed in terms of the commutator of the observables 
Q and B"{t), as 0bq(s) = {[Q,B"{s)])/ih (the super- 
script H denoting the Heisenberg picture with respect to 
the unperturbed dynamics.) Moreover, Kubo derived the 
general relation 



{QB''{t))^{B''it-ihf5)Q) 



(5) 



between differently ordered thermal correlation functions 
and deduced from it the celebrated quantum flu ctuation- 
dissipation theorem ( [Callen and Welton) 1951), reading: 



5fgg(a;) = {h/2i) coth{l3huj/2) 4>BQi(^) 



(6) 



where ^bq{^) — e*'^*5'Bg(s)ds, denotes the Fourier 
transform of the symmetrized, stationary equilibrium 
correlation function «'bq(s) = {QB" {s) + B"{s)Q)/2, 
and (Pbq{^) — /.!^ e*"^0BQ(s)ds the Fourier trans- 
form of the response function (pBQis)- Note that the 



Hanggi and Ingold ( 2005 i 



3 



fluctuation-dissipation theorem is valid also for many- 
particle systems independent of the respective particle 
statistics. Besides offering a unified and rigorous pic- 
ture of the fluctuation-dissipation theorem, the theory of 
Kubo also included other important advancements in the 
field of non-equilibrium thermodynamics. Specifically, 
we note the celebrated Onsager-Casimir reciprocity rela- 
tions ( jOnsager 1931a|b Casimir 1945). These relations 
state that, as a consequence of microreversibility, the ma- 
trix of transport coefficients that connects applied forces 
to so-called fluxes in a system close to equilibrium con- 
sists of a symmetric and an anti-symmetric block. The 
symmetric block couples forces and fluxes that have same 
parity under time-reversal and the antisymmetric block 
couples forces and fluxes that have different parity. 



Most importantly, the analysis of Kubo (1957) opened 



the possibility for a systematic advancement of response 
theory, allowing in particular to investigate the exis- 
tence of higher order fluctuation-dissipation relations, be- 
yond linear regime. This task was soon undertaken by 
Bernard and Callen (1959), who pointed out a hierar- 



chy of irreversible thermodynamic relationships. These 
higher order fluctuation dissipation relations were inves- 
tigated in detail by Stratonovich for Markovian system, 
and later by Efremov ( 1969 ) for non-Markovian systems, 
see ( Stratonovich 1992 Ch. I) and references therein. 



Even for arbitrary systems far from equilibrium the 
linear response to an applied force can likewise be re- 
lated to tailored two-point correlation functions of cor- 
responding stationary nonequilibrium fluctuations of the 
unde rlying unperturbed, stationary nonequilibri um sys- 
tem (jHanggT, T^TS", 'HanggTa^^ These 



authors coined the expression "fluctuation theorems" for 
these relations. As in the near thermal equilibrium case, 
also in this case higher order nonlinear response can be 
linked to corresponding higher order correlation functions 



of those nonequilibrium fluctuations ( Hanggi 1978 Prost 
et aLl 120091). 



At the same time, in the late seventies of the last 



century Bochkov and Kuzovlev (1977) provided a single 



compact classical expression that contains fluctuation re- 
lations of all orders for systems that are at thermal equi- 
librium when unperturbed. This expression, see Eq. ( 14) 



below, can be seen as a fully nonlinear, exact and uni- 
versal fluctuation relation. This IBochkov and KuzovlevI 
formula, Eq. ( |14[ ) below, soon turned out useful in ad- 
dressing the problem of connecting the deterministic and 
the stochastic descriptions of nonlinear dissipative sys- 



tems (Bochkov and Kuzovlev 1978 Hanggi 1982). 



As it often happens in physics, the most elegant, com- 
pact and universal relations, are consequences of general 
physical sym metries. In the case of [Bochkov and Ku-| 
zovlev (1977) the fluctuation relation follows from the 



time reversal invariance of the equations of microscopic 
motion, combined with the assumption that the system 
initially resides in thermal equilibrium described by the 
classical analogue of the Gibbs stat e, Eq. [Bochkov 



and Kuzovlev] ( |1977[ |1979[ |1981a|b[ ) proved Eq. ([l4|) be 



low for classical systems. Their derivation will be re- 
viewed in the next section. The quantum version, see Eq. 



( 55 ), was not r eported u ntil ve ry recently ( Andrieux and 
Gaspardl [20081). In Sec. ITTlC we shall discuss the funda- 



mental obstacles that pre vented , Bochkov and Kuzovlev 
( |1977[ |1979[ |1981a|b[ ) and |Stratonovich| ( |1994[ ) who also 
studied this very quantum problem. 

A new wave of activity in fluctuation relations was ini 
tiated by the works of Evans et al. ( 1993 ) and Gallavotti 



and Cohen (1995) on the statistics of the entropy pro- 



duced in non-equilibrium steady states, and of |Jarzynski| 
( 1997 ) on the statistics of work performed by a tran- 



sient, time-dependent perturbation. Since then, the field 
has generated grand interest and flourished considerably. 
The existing reviews on this topic mostly cover classi- 
cal fluctuat ion relations ([Jarzynski 2008 2011 Marconi 



eFoI] [20081 [Rondoni and Meji'a-MonasterioH2007[|Seifert 
2008'), while the comprehensive review by lEsposito et al. 



(j2009j) provides a solid, though in parts technical account 
of the state of the art of quantum fluctuation theorems. 
With this work we want to present a widely accessible 
introduction to quantum fluctuation relations, covering 
as well the most recent decisive advancements. Particu- 
larly, our emphasis will be on (i) their connection to the 
linear and non-linear response theory. Sec. [ll| (ii) the 
clarification of fundamental issues that relate to the no- 
tion of "work" , Sec. Ill (iii) the derivation of quantum 
fluctuation relations for both, closed and open quantum 
systems, Sec. [IV[ and [v| and also (iv) their impact for 



experimental applications and validation, Sec. VI 



II. NONLINEAR RESPONSE THEORY AND CLASSICAL 
FLUCTUATION RELATIONS 

A. Microreversibility of non autonomous classical systems 

Two ingredients are at the heart of fluctuation rela- 
tions. The flrst one concerns the initial condition of the 
system under study. This is supposed to be in thermal 
equilibrium described by a canonical distribution of the 
form of Eq. ([2|. It hence is of statistical nature. Its 
use and properties are discussed in many textbooks on 
statistical mechanics. The other ingredient, concerning 
the dynamics of the system is the principle of microre- 
versibility. This point needs some clarification since mi- 
croreversibility is customarily understood as a property 
of autonomous (i.e., non-driven) systems described by a 
time-independent Hamiltonian (Messiah 1962 Vol. 2, 
Ch. XV). On the contrary, here we are concerned with 
non-autonomous systems, governed by explicitly time- 
dependent Hamiltonians. In the following we will an- 
alyze this principle for classical systems in a way that 
at first glance may appear rather formal but will prove 
indispensable later on. The analogous discussion of the 
quantum case will be given next in Sec. |IV.A[ 

We deal here with a classical system characterized by a 
Hamiltonian that consists of an unperturbed part Hq (z) 



4 



and a perturbation —XtQ{z) due to an external force At 
that couples to the conjugate coordinate Q(z). Then the 
total system Hamiltonian becomes^ 



H{z,Xt)^Hoiz)-XtQiz), 



(7) 



where z = (q, p) denotes a point in the phase space of 
the considered system. In the following we assume that 
the force acts within a temporal interval set by a starting 
time and a final time r. The instantaneous force values 
At are specified by a function A, which we will refer to 
as the force protocol. In the sequel, it will turn out nec- 
essary to clearly distinguish between the function A and 
the value At that it takes at a particular instant of time 
t. 

For these systems the principle of microreversibility 
holds in the following sense. The solution of Hamilton's 
equations of motion assigns to each initial point in phase 
space Zq = (qo, Po) a point Zt at the later time t £ [0, r], 
which is specified by the values of the force in the order of 
their appearance within the considered time span. Hence, 
the position 



Zt = (pt,o[zo; A] 



(8) 



at time t is determined by the flow ipt,o[zo; A] which is a 
function of the initial point zq and a functional of the 
force protocol A."^ In a computer simulation one can 
invert the direction of time and let run the trajectory 
backwards without problem. Although, as experience ev- 
idences, it is impossible to actively revert the direction of 
time in any experiment, there is yet a way to run a time 
reversed trajectory in real time. For the sake of simplic- 
ity we assume that the Hamiltonian Hq is time reversal 
invariant, i.e. that it remains unchanged if the signs of 
momenta are reverted. Moreover we restrict ourselves to 
conjugate coordinates Q(z) that transform under time 
reversal with a definite parity eg 



Stratonovich 



= ±1 ^ 

(1994, Sec. 1.2. 3) showed that then the flow under the 
backward protocol A, with 



At — \t- 



(9) 



is related to the flow under the forward protocol A, via 
the relation 

<^t,o[zo;A] = e<^r-t,o[£Zr;eQA], (10) 
where e maps any phase space point z on its time reversed 



image ez = e(q, p) = (q, — p)- Equation (10) expresses 
the principle of microreversibility in driven systems. Its 
meaning is illustrated in Fig. [2] Particularly, it states 
that in order to trace back a trajectory, one has to reverse 
the sign of the velocity, as well as the temporal succession 
of the force values A, and the sign of force, A, if eq = — 1. 



^ The generalization to the case of several forces coupUng via dif- 
ferent conjugate coordinates is straightforward. 

^ Due to causality <pt,o[^0, ^] may of course only depend on the 
part of the protocol including times from t = up to time t. 



P 



£Zo 



<^t,o[zo; A] 



q 



¥'T-i,o[£Zr;£(3A] 



FIG. 2 Microreversibility for non-autonomous classical 
(Hamiltonian) systems. The initial condition zo evolves, un- 
der the protocol A, from time t — until time t to Zt = 
iy3f,o[zo; A] and until time t = r to Zt. The time-reversed final 
condition ezr evolves, under the protocol eg A from time i = 
until T — t to <p-r-t,o[£z-r;eQX] = s<pt,o[zo; X], Eq. (10 1, and 



until time t — t to the time-reversed initial condition ezq. 



B. Bochkov-Kuzovlev approach 

We consider a phase-space function B{z) which has 
a definite parity under time reversal eb = ±1, i-e. 
B{ez) — eBB{z). Let Bt = S((^t,o[zo; A]) denote its 
temporal evolution. Depending on the initial condition 
Zq different trajectories Bt are realized. Under the above 
stated assumption that at time t = the system is pre- 
pared in a Gibbs equilibrium, the initial conditions are 
randomly sampled from the distribution 



po(zo) - e-^^"(^°)/Zo, 



(11) 



with Zq — J dzQe~'^^'>''^°\ Consequently the trajectory 
Bt becomes a random quantity. Next we introduce the 
quantity: 



Wo[zo;A] 



dtXtQt, 



(12) 



where Qt is the time derivative of Qt — Q{i pt,o\zQ; A]). 



From Hamilton's equations it follows that (Jarzynski 
20071): 



Wo[zo; A] = iJo(</5T,o[zo; A]) - i?o(zo)- 



(13) 



Therefore, we interpret Wq as the work injected in the 
system described by Hq during the action of the force 



5 



protocol.^ The central finding of Bochkov and Kuzovlev 



( 1977 ) is a formal relation between the generating func- 



tional for multi-time correlation functions of the phase 
space functions Bt and Qt and the generating functional 
for the time-reversed multi-time auto-correlation func- 
tions of Bt , reading 



releasing events which, in fact, must happen with cer- 
tainty in order that Eq. (151 holds if the average work 
is larger than zero. Equation ( 16 1 may be regarded as 



a microscopic manifestation of the Second Law of ther- 
modynamics. For this reason .Stratonovich ( 1994 Sec. 



(14) 



where u-r is an arbitrary test-function, ut = Ur~t is its 
temporal reverse, and the average denoted by (•) is taken 



with respect to the Gibbs distribution po of Eq. (11). On 
the left hand side. 



are 



the time evolutions of Bt and Qt 
governed by the full Hamiltonian ([t]) in presence of the 
forward protocol as indicated by the subscript A while on 
the right hand side the dynamics is determined by the 
time-reversed protocol indicated by the subscript £qA. 
The derivation of Eq. (14 1, which is based on the mi- 



croreversibility principle, Eq. (10), is given in Appendix 
[A| The importance of Eq. ( 14 ) lies in the fact that it 
contains the Onsager reciprocity relations and fluctua- 
tion relations of all orders within a single compact for- 
mula (Bochkov and Kuzovlev 1977). These relations 



may be obtained by means of functional derivatives of 
both sides of the Eq. ( 14 ), of various orders, with respect 



1.2.4) referred to it as the H-theorem. We recapitulate 
that only two ingredients - initial Gibbsian equilibrium 
and microreversibility of the dynamics - have led to Eq. 
( [l4| ). In conclusion, this relation not only contains linear 
and nonlinear response theory, but also the second law 
of thermodynamics. 

The complete information about the statistics is con- 
tained in the work probability density function (pdf) 
Po[W^o;A]. The only random element entering the work, 
Eq. ( 13 ), is the initial phase point Zq which is distributed 
according to Eq. (11). Therefore po ; A] may formally 
be expressed as 

Po[W^o;A] = y"dzopo(zo)<5[M^o-i^o(zr)+i^o(zo)], (17) 

where 5 denotes Dirac's delta function. The hmctional 
dependence of po[W^o;A] on A is contained in the term 
•Zr = (pr,o[zo;^]- Using the microreversibility principle, 
Eq. (10), one obtains the following fluctuation relation 



to the force-field A and the test-field u at vanishing fields 
A = M = 0. The classical limit of the I C alien and Wel-l 
(1951) fluctuation-dissipation theorem, Eq. ([6]), for 



ton 



Po[Wo;A]^ 

Po[-Wo]£q\] 



(18) 



instance, is obtained by differentiation with respect to m, in a way analogous to the derivation of Eq. ([14|). We 



followed by a differentiation with respect to A (Bochkov 



and Kuzovlev 1977), both at vanishing fields u and A. 



Another remarkable identity is achieved from Eq. ( 14 ) 
by putting u = 0, but leaving the force A finite 
yields the Bochkov-Kuzovlev equality, reading 



This 



shall refer to this relation as the Bochkov-Kuzovlev work 
fluctuation relation, although it was not explicitly given 
by Bochkov and Kuzovlev, but was only recently ob- 
tained by Horowitz and Jarzynski ( 2007 ) . This equation 



1. 



(15) 



has a profound physical meaning. Consider a positive 
work Wo > 0. Then Eq. (18) says that the probabil- 



In other words, for any system that initially stays in ther- 
mal equilibrium at a temperature T = 1/ [ksP) the work. 



Eq. ( 12 ) , done on the system by an external force is a ran- 



dom quantity with an "exponential expectation value" 
{e~^'^'')\ that is independent of any detail of the system 
and the force acting on it. This of course does not hold 
for the individual moments of work. Since the exponen- 
tial function is concave, a direct consequence of Eq. ( 15 ) 
is 



ity that this work is injected into the system is larger 
by the factor e^'^° than the probability that the same 
work is absorbed under the reversed forcing: energy con- 
suming processes are exponentially more probable than 
energy releasing processes. Thus, Eq. ( 18 ) expresses the 



second law of thermodynamics at a very detailed level 
which quantifles the relative frequency of energy releas- 
ing processes. By multiplying both sides of Eq. ( 18 ) by 
Po[—WQ;eQ\]e~^^° and integrating over VFq, one recov- 
ers the Bochkov-Kuzovlev identity, Eq. (15). 



(Wo) A > 0. 



(16) 



That is, on average, a driven Hamiltonian system may 
only absorb energy if it is perturbed out of thermal equi- 
librium. This does not exclude the existence of energy 



Following || jarzynski| |2007| we refer to Wq as the exclusive 
work, to distinguish from the inclusive work W = H{zr, A^) — 
-ff(zo,Ao), Eq. ( |19[ ), which accounts also for the coupling be- 
tween the external source and the system. Wc will come back 
later to these two definitions of work in Sec. IIII.AI 



C. Jarzynski approach 

An alternative deflnition of work is based on the com- 
parison of the total Hamiltonians at the end and the be- 
ginning of a force protocol, leading to the notion of "in- 
clusive work" in contrast to the "exclusive work" defined 



in Eq. ( 13 ). The latter equals the energy difference refer- 
ring to the unperturbed Hamiltonian Hq. Accordingly, 
the inclusive work is the difference of the total Hamilto- 
nians at the final time t = t and the initial time t = Q: 



W[zo;X]=H{zr,Xr)-H{zoAo) 



(19) 



6 



In terms of the force At and the conjugate coordinate Qt, 
the inclusive work is expressed as'"^: 



W^[zo; A] = / dtX 



(20) 



dH{zt,Xt) 
• dXt 

dtXtQt 
Wo[zo;X]-XtQt\l 



For the sake of simpUcity we confine ourselves to the case 
of an even conjugate coordinate Q. In the corresponding 
way, as described in AppendixjXj we obtain the following 
relation between generating functionals of forward and 
backward processes in analogy to Eq. (14 1, reading 



ZjXr) 
ZiXo) 



/q dtuteBBt\^_ 



[e-'o 



(21) 



While on the left hand side the time evolution is con- 
trolled by the forward protocol A and the average is per- 
formed with respect to the initial thermal distribution 
/9^(z, Ao), on the right hand side the time evolution is gov- 
erned by the reversed protocol A and averaged over the 
reference equilibrium state p^(z,Ao) 



— p^{z, Xr)- Here 



p^(z,AO = e-^^(-^*)/Z(A,) 



(22) 



formally describes thermal equilibrium of a system with 
the Hamiltonian H{z, At) at the inverse temperature /3. 
The partition function Z(Xt) is defined accordingly as 
Z{Xt) = J dze~^"^^'^'\ Note that in general the ref- 
erence state p^(z,At) is different from the actual phase 
space distribution reached under the action of the pro- 
tocol A at time t, i.e., p{z,t) — p''((^^g [z; A], Aq), where 
(fi^ q[z;X] denotes the point in phase space that evolves 
to z in the time to t under the action of A. 
Setting u = we obtain 



where 



AF = F(A^) - F{Xo) = -/31n 



ZjXr) 

Z{Xo)- 



(23) 



(24) 



is the free energy difference between the reference state 
p^{z, At) and the initial equilibrium state p^{zXo)- As a 
consequence of Eq. (p3| we have 



(25) 



{W)x > AF, 



which is yet another expression of the second law of ther- 
modynamics. Equation (23) was first put forward by 



Jarzynski ( 1997 ) , and is commonly referred to in the lit- 



erature as the "Jarzynski equality" 



For a further discussion of inclusive and exclusive work we refer 

to Sect. \nrK\ 



In close analogy to the Bochkov-Kuzovlev approach 
the pdf of the inclusive work can be formally expressed 
as 

p[W;X]= J dzopizo,Xo)6[W - Hizr,Xr) + Hizo,Xo)]. 

(26) 

Its Fourier transform defines the characteristic function 
of work: 

G[u; A] = / dWe'^'^plW; A] 



dzoe^"^^^^^'^^^^"'-^"'^"''^e^^^^^°'^°^ /Z{Xo) 



= / dzo exp 



iu I dtXt 





dHjzt^Xt) 
dXt 



Z{Xo) ■ 



(27) 



Using the microreversibility principle, Eq. (10), we 



obtain in a way similar to Eq. ( 18 ) the (inclusive) work 
fiuctuation relation: 



P[W;X]^ 
p[-W;X] 



(W-AF) 



(28) 



where the probability p[—W; X] refers to the backward 
process which for the inclusive work has to be determined 
with reference to the initial thermal state p^{z, Xr)- First 
put forward by Crooks ( 1999 ) , Eq. ( 28 ) is commonly re- 



ferred to in literature as the "Crooks fluctuation theo- 
rem". The Jarzynski equality, Eq. (23), is obtained by 

-W;X]e-^^ and 



23 28) continue to 



multiplying both sides of Eq. ( 28 ) by p[ 
integrating over W. Equations (21 
hold also when Q is odd under time reversal, with the 
provision that A is replaced by —A. 

We here point out the salient fact that, within the in- 
clusive approach, a connection is established between the 
nonequilibrium work W and the difference of free ener- 
gies AF, of the corresponding equilibri um s tates p^{z, Xt) 
and p^{z,Xo). Most remarkably, Eq. (25) says that the 



average (inclusive) work is always larger than or equal 
to the free energy difference, no matter the form of the 
protocol A; even more surprising is the content of Eq. 
( 23 ) saying that the equilibrium free energy difference 



may be inferred by measurements of nonequilibrium work 
in many realizations of the forcing experiment (Jarzyn- 
ski[ 1997). This is similar in spirit to the fluctuation- 
dissipation theorem, also connecting an equilibrium prop- 
erty (the fluctuations) , to a non-equilibrium one (the lin- 
ear response), with the major difference that Eq. (23 1 is 



an exact result, whereas the fluctuation-dissipation the- 
orem holds only to flrst order in the perturbation. Note 
that as a consequence of Eq. ( 28 ) the forward and back- 



ward pdf's of exclusive work take on the same value at 
W = AF. This property has been used in experiments 



( Collin et al. 


2005 


Douarche et al. 


2005 iLiphardt et al. 


2002p in order to determine free energy differences from 



nonequilibrium measurements of work. Equations (23 
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281 have further been employed to develop efficient nu- 



merical methods for the estimation of free energies ( 


Hahn 


and Then, 2009 2010 Jarzynski| |2002| Minh and Adib| 


2008 


Vaikuntanathan and Jarzynski 


2008 


1. 



Both the Crooks fluctuation theorem, Eq. (28), and 



the Jarzynski equality, Eq. (23), continue to hold for any 
time dependent Haniiltonian -ff (z. At) without restriction 
to Hamiltonians of the form in Eq. ([?]). Indeed no re- 
striction of the form in Eg. ([7 ) was imposed in the sem- 
inal paper by Jarzynski (1997). In the original works of 
Jarzynskil (11997^ and ICrooksl ( 19991), Eqs. (pSl) and (IM 



were obtained directly, without passing through the more 
general formula in Eq. (21 ). Notably, neither these sem- 



inal papers, nor the subsequent literature, refer to such 
general functional identities as Eq. (21). We introduced 



them here to emphasize the connection between the re- 
cent results, Eqs. (23) and (28), with the older results 
of iBochkov and Kuzovlevl (119771), Eqs. (flSl pi). The 



latter ones were practically ignored, or sometimes misin- 
terpreted as special instances of the former ones for the 
case of cyclic protocols {AF = 0), by those working in 
the field of non-equilibrium work fluctuations. Only re- 
cently I Jarzynski I ([2007j) pointed out the differences and 
analogies between the inclusive and exclusive approaches. 



III. FUNDAMENTAL ISSUES 

A. Inclusive, exclusive and dissipated work 

As we evidenced in the previo us se ction, the s tudies 
of Bochkov and Kuzovlev (1977) and Jarzynski (1997) 



are based on different definitions of work, Eqs . (|13 



19), reflecting two different viewpoints (Jarzynski 2007). 



From the "exclusive" viewpoint of |Bochkov and Kuzovlev] 
( 1977 ) the change in the energy Hq of the unforced sys- 
tem is considered, thus the forcing term {—XtQ) of the 
total Hamiltonian is not included in the computation of 
work. From the "inclusive" point of view the definition 
of work, Eq. (19), is based on the change of the to- 



tal energy H including the forcing term {—XtQ)- In ex- 
periments and practical applications of fluctuation rela- 
tions, special care must be paid in properly identifying 
the measured work with either the inclusive (W) or ex- 
clusive (Wq) work, bearing in mind that A represents the 
prescribed parameter progression and Q is the measured 
conjugate coordinate. 



The experiment of Douarche et at. ( 2005 ) is very well 



suited to illustrate this point. In that experiment a pre- 
scribed torque Mt was applied to a torsion pendulum 
whose angular displacement 9t was continuously moni- 
tored. The Hamiltonian of the system is 



H{y,pe,e,Mt) =i?B(y) + i?ss(y,P9 

-f — H 

21 2 



Mt9, (29) 



where pg is the canonical momentum conjugate to 9, 
Hb (y) is the Hamiltonian of the thermal bath to which 



the pendulum is coupled via the Hamiltonian Hsbi and 
y is a point in the bath phase space. Using the defini- 
tions of inclusive and exclusive work, Eqs. (|T2l[20|, and 



, Eqs. ([12] 

noticing that M plays the role of A and 9 that of Q, we 
find in this case W = - J 9Mdt and IVq = / M9At. 

Note that the work W = - ^ 9Mdt, obtained by mon- 
itoring the pendulum degree of freedom only, amounts to 
the energy change of the total pendulum-|-batli system. 
This is true in general (Jarzynski 2004). Writing the 
total Hamiltonian as 

i/(x,y,At) = i/5(x,A0+ffBs(x,y) + HB(y), (30) 

with Hs{x,Xt) being the Hamiltonian of the system of 
interest, one obtains 



9Xt 



, dH 
dt 



J TT 

at-^w, (31) 



because Hbs and Hb do not depend on time, and as a 
consequence of Hamilton's equations of motion dH/dt = 
dH/dt. 

Introducing the notation M^diss — W — AF, for the 
dissipated work, one deduces that the Jarzynski equality 
can be re-expressed in a way that looks exactly as the 
Bochkov-Kuzovlev identity, namely: 



1. 



(32) 



This might let one believe that the dissipated work coin- 
cides with Wq. This, however, would be incorrect. As dis- 



cussed in (Jarzynski 2007) and explicitly demonstrated 
by Campisi et al. (2011a) Wq and Wdiss constitute dis- 
tinct stochastic quantities with different pdf's. The in- 
clusive, exclusive and dissipated work coincide only in the 
case of cyclic forcing X^ = Aq (Campisi et al. 2011a). 



B. The problem of gauge freedom 

We point out that the inclusive work W, and free en- 
ergy difference AF, as defined in Eqs. (19, 24), are - 
to use the expression coined by [Cohen- Tannoudji et al. 
(1977) - not "true physical quantities." That is to say 



they are not invariant under gauge transformations that 
lead to a time-dependent shift of the energy reference 
point. To elucidate this, consider a mechanical sys- 
tem whose dynamics are governed by the Hamiltonian 
H{z,Xt). The new Hamiltonian 



H'{z,Xt)^H{z,Xt)+g{Xt 



(33) 



where g{Xt) is an arbitrary function of the time depen- 
dent force, generates the same equations of motion as H. 
However, the work W = H'{zT,Xr) — H'{zo,Xo) that 
one obtains from this Hamiltonian differs from the one 
that follows from 77, Eq. ([l9|: W' ^ W + g{Xr) - g{Xo). 
Likewise we have, for the free energy difference AF' — 
AF + g{Xr) — giXo) . Evidently the Jarzynski equality, Eq. 
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(231, is invariant under such gauge transformations, be- C. Work is not a quantum observable 



cause the term g(Ar) — <7(Ao) appearing on both sides of 
the identity in the primed gauge, would cancel; explicitly 
this reads: 



(34) 



Thus, there is no fundamental problem associated with 
the gauge freedom. 

However one must be aware that, in each particular 
experiment, the very way by which the work is mea- 
sured implies a specific gauge. Consider for example 



the torsion pendulum experiment of ( jPouarche et al. 
|2005p . The inclusive work was computed as: W = 
— J 6Mdt. The condition that this measured work is 
related to the Hamiltonian of Eq. ( p9| via the rela- 
tion W — i?(z7-,Ar) — H{zq,Xo) , Eq. p^, is equiv- 



alent to -J^0Mdt = J^{dH/dM)Mdt, see Eq. (jsi]) 
If this is required for all r then the stricter condition 
dH/dM — ^9 is implied, restricting the remaining gauge 
freedom to the choice of a constant function g. This resid- 
ual freedom however is not important as it does neither 
affect work nor free energy. We now consider a different 
experimental setup where the support to which the pen- 
dulum is attached is rotated in a prescribed way accord- 
ing to a protocol a^, specifying the angular position of 
the support with respect to the lab frame. The dynamics 
of the pendulum are now described by the Hamiltonian 



Iuj'ate + g{at). (35) 



If the work W — J Nadt done by the elastic torque 
N — /w^(a — 9) on the support is recorded then the 
requirement dH/da = N singles out the gauge g{af) = 
+ const, leaving only the freedom to chose the 
unimportant constant. Note that when Mt = luPat, the 
pendulum obeys exactly the same equations of motion 
in the two examples above, Eqs. (29 35 1. The gauge is 



irrelevant for the law of motion but is essential for the 
energy-Hamiltonian connection.^ 

The issue of gauge freedom was first pointed out by 
Vilar and Rubi (2008c), who questioned whether a con- 



nection between work and Hamiltonian may actually ex- 
ist. Since then this topic had been highly debated,'' but 
neither the gauge invariance of fluctuation relations nor 
the fact that different experimental setups imply different 
gauges were clearly recognized before. 



bee also |Kobe| ( |l98lt , in the context of non-relativistic electro- 
dynamics. 

See |Adib| |2009||; |Chen l |)2008alb| l2009|l; |Crooks| ||2009|l; 
[Horowitz and Jar2yns'kilT 2008); Pcliti ( 2008a'^ n'Vilar and Rubi] 
( |2008a|b|c[ |; |Zimanyi and Silbey, ( |2009) ). 



Thus far we have reviewed the general approach to 
work fluctuation relations for classical systems. The 
question then naturally arises of how to treat the quan- 
tum case. Obviously, the Hamilton function H{z,Xt) is 
to be replaced by the Hamilton operator 'H(At), Eq. (|3|. 
The probability density p{z, At) is then replaced by the 
density matrix g{Xt), reading 



/2{Xt), 



(36) 



where Z{Xt) = Tre~'^^^^'-* is the partition function and 
Tr denotes the trace over the system Hilbert space. The 
free energy is obtained from the partition function in 
the same way as for classical systems, i.e., F{Xt) — 
— InZ(At). Less obvious is the deflnition of work in 
quantum mechanics. Originally, [Bochkov and Kuzovlev| 
( 1977 ) defined the exclusive quantum work, in analogy 



with the classical expression, Eqs. (12 13 1, as the opera- 
tor Wo = Jq dtXt = — Ho , where the superscript 
H denotes the Heisenberg picture: 



^ul„[X]BUtM. 



(37) 



Here B is an operator in the Schrodinger picture and 
?7(^o[A] is the unitary time evolution operator governed 
by the Schrodinger equation 

,n^d^^U{Xt)UM C^o,o[A]=l, (38) 
ot 

with 1 denoting the identity operator. We use the no- 
tation J7t,o[A] to emphasize that, like the classical evo- 
lution (/3t_o[z;A] of Eq. ([8]), the quantum evolution op- 
erator is a functional of the protocol A.^ The time 
derivative Qf is determined by the Heisenberg equa- 
tion. In case of a time-independent operator Q it be- 
com es gf = ^[?^f(At),gf| A^ 
,Bochkov and Kuzovlev] ( |1977 ) were not able to pro- 
vide any quantum analog of their fluctuation relations, 
Eqs. (14 151, with the classical work Wq replaced by the 



operator >V( 



Yukawa (20001 and Allahverdyan and Nieuwenhuizen 



( 2005 ) arrived at a similar conclusion when attempting 
to deflne an inclusive work operator by W — V.^ {X^) — 
■H(Ao). According to this definition the exponentiated 
work (e-^^) = Trgoe-P^ = (e-'3[«?(^-)-'^(^o)]) agrees 
with e~^^^ only if 'H{Xt) commutes at different times 
['H{Xt),7i{Xr)] = for any t,T. This could lead to the 
premature conclusion that there exists no direct quantum 
analog of the Bochkov-Kuzovlev and the Jarzynski iden- 



tities, Eqs. (15 231, (Allahverdyan and Nieuwenhuizen 
20051). 



Due to causality (7t,o[-^] of course only depend on the part 

of the protocol including times from up to t. 
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Based on the works by Kurchan (2000) and Tasaki 



(2000), [Talkner et al. (2007) and [Talkner and Hanggi 



(2007) demonstrated that this conclusion is based on an 



erroneous assumption. They pointed out that work char- 
acterizes a process, rather than a state of the system; this 
is also an obvious observation from thermodynamics: un- 
like internal energy, work is not a state function (its dif- 
ferential is not exact). Consequently, work cannot be 
represented by a Hermitean operator whose eigenvalues 
can be determined in a single projective measurement. 
In contrast, the energy 'H{Xt) (or Hq, when the exclusive 
viewpoint is adopted) must be measured twice first at the 
initial time t = and again at the final time time t — t. 

The difference of the outcomes of these two measure- 
ments then yields the work performed on th e system in 
a particular realization ( Talkner et al. 2007 ) . That is, if 
at time t = the eigenvalue E:^° of H(Ao) and later, at 
t = T, the eigenvalue of T-L{\t) were obtained,^ the 
measured (inclusive) work becomes: 



A. Microreversibility of non autonomous quantum systems 

The principle of microreversibility is introduced and 
discussed in quantum mechanics textbooks for au- 
tonomous (i.e., non-driven) quantum systems (Messiah 
1962). As in the classical case, however, this principle 



continues to hold in a more general sense also for non- 
autonomous quantum systems. In this case it can be 
expressed as: 



(40) 



wher e 6 is the quan tum mechanical time reversal oper- 



ator (Messiah 1962).^ Note that the presence of the 



protocol A and its time reversed image A distinguishes 
this generalized version from the standard form of mi- 
croreversibility for autonomous systems. The principle 
of microreversibility, Eq. ( 40 ) , holds under the assump- 



tion that at any time t the Hamiltonian is invariant under 
time reversal,^"'^ that is: 



Et 



(39) 



Equation ( 39 ) represents the quantum version of the 
classical inclusive work, Eq. ( 19 ). In contrast to the clas- 



sical case, this energy difference, which yields the work 
performed in a single realization of the protocol, cannot 
be expressed in the form of an integrated power, as in 
Eq. ([20l). 



The quantum version of the exclusive work, Eq. 
is wq = e„ 



(131, 



2011a), where now e/ are 



- e„ ( Campisi et al. 
the eigenvalues of T-Lq. As we will demonstrate in the next 
section, with these definitions of work straightforward 
quantum analogs of the Bochkov-Kuzovlev results, Eqs. 
( 14 15 18 ) and of their inclusive viewpoint counterparts, 
Eqs. (1211 1231 [281) can be derived. 



IV. QUANTUM WORK FLUCTUATION RELATIONS 

Armed with all the proper mathematical definitions of 
nonequilibrium quantum mechanical work, Eq. (39), we 



next confidently embark on the study of work fluctua- 
tion relations in quantum systems. As in the classical 
case, also in the quantum case one needs to be careful 
in properly identifying the exclusive and inclusive work, 
and must be aware of the gauge freedom issue. In the fol- 
lowing we shall adopt, except when otherwise explicitly 
stated, the inclusive viewpoint. The two fundamental 
ingredients for the development of the theory are, in the 
quantum case like in the classical case, the canonical form 
of equilibrium and microreversibility. 



(41) 



A derivation of Eq. ( 40 ) is presented in Appendix p3l 



See also ( Andrieux and Gaspard 2008 ) for an alternative 



derivation. 

In order to better understand the physics behind Eq 
(40) we rewrite it as: ?7t,o[A] = Q^Ur 
where we used the concatenation 



*,o[A]ei7,.o[A], 
rule C/t.r[A] = 



C^i,o[A]?7o,r [A], and the inverse C/o,r[A] = t/Tro[A] of the 
propagator Ut^X]. Applying it to a pure state \i), and 
multiplying by Q from the left, we obtain: 



(42) 



where we introduced the notations IV^t) = ^^t,o[A]|j) and 
I/) = [/r,o[A]|«)- Equation (42) says that, under the evo- 



lution generated by the reversed protocol A the time re- 
versed final state, 6|/), evolves between time and r — i, 
to Q\'ipt)- This is illustrated in Fig. [3j As for the classical 
case, in order to trace a non autonomous system back to 
its initial state, one needs not only to invert the momenta 
(applying Q), but also to invert the temporal sequence 
of Hamiltonian values. 



B. The work probability density function 

We consider a system described by the Hamiltonian 
'H(At) initially prepared in the canonical state Q{\a) = 



For a formal definition of these eigenvalues see Eq. J43| below. 



Under the action of coordinates transform evenly, whereas 

linear and angular momenta, as well as spins change sign. In the 

coordinate representation, in absence of spin degrees of freedom, 

the operator © is the complex conjugation: Q>ip = ip* . 

In presence of external magnetic fields the direction of these fields 

has also to be inverted in the same way as in the autonomous 

case. 
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\f) 






ei/) 









FIG. 3 (Color online) Microreversibility for non autonomous 
quantum systems. The normalized initial condition \i) 
evolves, under the unitary time evolution generated by H(At), 
from time t = until t to \tpt} ~ Ut,o[X]\i) and until time t = t 
to I/). The time-reversed final condition &\f) evolves, under 
the unitary evolution generated by A f rom time t — until 
T-t to (7r-t,o[A]01/) = 0\i>t), Eq. (421, and until time t = r 
to the time-reversed initial condition 0\i) . The motion occurs 
onto the hypersphere of unitary radius in the Hilbert space. 



e '^'^ /Z{Xo). The instantaneous eigenvalues of 'H{Xt) 
are denoted by E^* , and the corresponding instantaneous 



eigenstates by |'0^*^) 



(43) 



The symbol n labels the quantum number specifying 
the energy eigenvalues and 7 denotes all further quan- 
tum numbers, necessary to specify an energy eigen- 
state in case of g„-fold degeneracy. We emphasize that 
the instantaneous eigenvalue equation ( 43 1 must not be 



confused with the Schrodinger equation ihdt\ipit)) = 
H{Xt)\4'{t)) ■ The instantaneous eigenfunctions result- 
ing from (|43| in particular are not solutions of the 
Schrodinger equation. 

At t = the first measurement of H{Xo) is performed, 
with outcome . This occurs with probability 



- g„e-^^"7Z(Ao) 



(44) 



According to the postulates of quantum mechanics, im- 
mediately after the measurement of the energy the 
system is found in the state: 



(45) 



= J2-y li'n^-y) ii^n'l'yl IS the projcctor onto the 



where ZI^" 

eigenspace spanned by the eigenvectors belonging to the 
eigenvalue E^° . The system is assumed to be thermally 



isolated at any time t > 0, so that its evolution is deter- 



mined by the unitary operator Ut.olX], Eq. (38), hence it 
evolves according to 



Quit) = Ut,o[^]g.nUUx]. 



(46) 



At time r a second measurement of HlXr) yielding the 
eigenvalue E^ with probability 



(47) 



is performed. The pdf to observe the work w is thus given 
by: 



p[w; A] = ^ Siw - [E^- - i?^])p™|„[AK. 



(48) 



The work pdf has been calculated explicitly for a forced 
harmonic oscillator ( Talkner a/. 2008a 2009a I and for 



a parametric oscillator with varying frequency (Deffner 
et al.\ 120101 iDeffner and Lutzl 120081). 



C. The characteristic function of work 

The characteristic function of work G[u;X] is defined 
as in the classical case as the Fourier transform of the 
work probability density function 



G[u; X] = / du;e™>[u;; A] 



(49) 



Like p[w;A], G[m;A] contains full information regarding 
the statistics of the random variable w. ITalkner et al.\ 



( 2007 ) showed that the work pdf has the form of a two- 



time non-stationary quantum correlation function; i.e., 

G[w;A] = (e^^^^^-^e-^^^^o)) (50) 
^ Tre™«?(^^)e-(^"+««(^°)/Z(Ao) , 



where the average symbol stands for quantum expecta- 
tion over the initial state density matrix g{Xo), Eq. (36 ), 
(B) = Trp(Ao)S, and the superscript H denotes 



Heisenberg picture, i.e., (Xr) — ulQ[X]'H{Xr)UT,o[X]. 



Equation ( 50 ) was derived first by Talkner et al. ( 2007 ) 
in the case of nondegener ate 'H(At) and later generalized 
by Talkner et al. (2008b I to the case of possibly degen- 



erate HiXt 

The product of the two exponential operators 
giu'H^ (A^)g-iu«(Ao) (jg^j^ combined into a single expo- 
nent under the protection of the time ordering operator 



The derivation in [Talkner et~al] l |2008b| l is more general in that 
it does not assume any special form of the initial state, thus 
allowing the study, e.g., of microcanonical fluctuation relations. 
The formal expression, Eq. ( |50[ l, remains valid for any initial 
state Q, with the provision that the average is taken with respect 
to § = Ifn" QTln° representing the diagonal part of p in the 
eigenbasis of 'H(Ao). 
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T to yield e™«?(^-)e-™^(^") = T e™[«?(^-)-«(^o)l . In and (|Stratonovichl|1994[) employed yet different ordering 



this way one may convert the characteristic function of 
work to a form that is analogous to the corresponding 
classical expression, Eq. (p7|), 



G{u-\\ = Trre"'^?(^-)-^(^°)e-^«(^°V2o (51) 



= TrTexp 



lu I dtXt- „^ 
oXt 



/2o 



The second line follows from the fact that the total time- 
derivative of the Hamiltonian in the Heisenberg picture 
coincides with its partial derivative. 

As a consequence of quantum microreversibility, Eq. 



(401, the characteristic function of work obeys the fol- 



lowing important symmetry relation (see Appendix O) 



Z{\q)G[u- A] = Z{\r)G[-u + ip- A] 



(52) 



By applying the inverse Fourier transform and using 
Z(At) = ^T-f,~m{M) = e-'3^(^*) one ends up with the 
quantum version of the Crooks fluctuation theorem in 
Eq. (pSl): 



p[-w; A] 



Ji(w-AF) 



(53) 



This result was first accomplished by Tasaki' ('2000') and 



Kurchan| ( |2000[ ). Later |Talkner and Hanggi, ( ,2007| ) gave a 
systematic derivation based on the characteristic function 
of work. The quantum Jarzynski equality 



-PAF 



(54) 



follows by multiplying both sides by p[— w; A]e-^'" and 
integrating over w. Given the fact that the characteristic 
function is determined by a two-time quantum correla- 
tion function rather than by a single time expectation 
value is another clear indication that work is not an ob- 
servable but instead characterizes a process. 

As discussed in (iCampisi et al. 2010a) the Tasaki- 
Crooks relation, 



Eq. 



(53) and the quantum version of 
the Jarzynski equality, Eq. (54), continue to hold even 



if further projective measurements of any observable A 
are performed within the protocol duration (0, r ). These 
measurements, however do alter the work pdf (Campisi 



et al. 2011b). 



D. Quantum generating functional 

The Jarzynski equality can also immediately been ob- 
tained from the characteristic function by setting u — if3, 
in Eq. (l50| (iTalkner et al.\ 120071). In order to obtain 



this result it is important that the Hamiltonian opera- 
tors at initial and final times enter into the characteristic 
function, Eq. ( 50 ) , as arguments of two factorizing ex- 
ponential functions, i.e. in the form e~ 



In general, this of course is different from a single ex- 
ponential e~^['^ (-'^t)-'H(Ao)]_ In the definitions of gener- 
ating functionals, Bochkov and Kuzovlev (1977 1981a 



prescriptions which though do not lead to the Jarzynski 
equality. In order to maintain the structure of the clas- 
sical generating functional, Eq. (21), also for quantum 



systems the classical exponentiated work e~^^ has to be 
replaced by the product of exponentials as it appears in 
the characteristic function of work, Eq. (50). This then 



leads to a desired generating functional relation 



exp 



dsutB^ 



g-/3«"(At)g/3H(Ao) 



exp 



dtutEisB^ 



-I3AF 



(55) 



where B is an observable with definite parity eg (i.e., 
QBQ^ = EbB), B^ denotes the observable B in the 
Heisenberg representation, Eq. (37), Ut is a real func- 
tion, and Ut — Ur-t- This can be proved by using the 
quantum microreversibility principle, Eq. (40), in a sim 
ilar way as in the classical derivation, Eq. (XT) 
The derivation of Eq 



( |55[ ) was provided by Andrieux 
land Gaspard (2008), who therefrom also recovered the 
formula of Kubo, Eq. ([4]), and the Onsager-Casimir reci- 
procity relations. ^'^ 

These relations are obtained by means of functional 
derivatives of Eq, 
At, and test fields, ut, at At — ut 



(55) with respect to the force fields, 
( Andrieux and Gas- 



pard 2008). Relations and symmetries for higher order 



response functions follow in an analogous way as in the 
classical case, Eq. (14), by means of higher order func- 
tional derivatives with respect to the force field Af. Such 
relations were investigated experimentally by [Nakamura 
et_al.\ (|2010[); [Nakamura et al.\ ( |2011[ ), see Sec. |V) Eq 



(89) 



Within the exclusive viewpoint approach, the coun- 
terparts of Eqs. (50 52l 53 54 55), are obtained by 



replacing, T-L^ with "H^ ^{S^ with "Ho, -Z(At) with 
Zq = Tre~^^°, and setting accordingly AF to zero 



(Campisi et al. 2011a 



E. l\/licroreversibility, conditional probabilities and entropy 

For a Hamiltonian 'H{Xt) with nondegenerate instan- 
taneous spectrum for all times t and instantaneous eigen- 
vectors IV'n') the conditional probability p„i|n['^]i Eq- 
(47), is given by the simple expression: Pmini-^] = 
I^T-,o['^]|'/'^'')P- As a consequence of the assumed 
invariance of the Hamiltonian, Eq. (41), the eigenstates 



Andrieux and Gaspard' ('2008') also allowed for a possible depen- 
dence of the Hamiltonian on a magnetic field "H = 'H(At,B). 
Then it is meant that the dynamical evolution of B in the right 
hand side is governed by the Hamiltonian T-Li^Xt, — B), i.e., besides 
inverting the protocol, the magnetic field needs to be inverted as 
well. 
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I tpn* ) ) are invariant under the action of the the time re- 
versal operator up to a phase, O|'0^*) = e**^" IV'n')- Then, 
expressing microreversibility as f/T-^oW = ©^t/o.rl^]©, see 
Eq. ( 40 ) , the following symmetry relation is obtained for 



the conditional probabilities: 

Pm I n 



(56) 



Note the exchanged position of m and n in the two sides 
of this equation. From Eq. (56) the Tasaki-Crooks fluc- 



tuation theorem is readily obtained for a canonical initial 
state, using Eq. (48 1. 



Using instead an initial microcanonical state at energy 
E, described by the density matrix^^ 

Po{E)^d[E-H{Xo)]/Lo{E,Xo), (57) 

where w^B, At) = Tr S{E - H(Xt)), we obtain ([Talkner 



et al. 2008b): 



P[E,W;X] ^ 
p[E + W, -W; X] 



^lS(E+W,Xr}-S(EM)]/kE 



(58) 



where S{E,Xt) — fc^ lnaj(i?, At), denotes Boltzmann's 
thermodynamic equilibrium entropy. The corresponding 
classical derivation was provided bylCleuren et al. ( 2006 1 . 



A classical microcanonical version of the Jarzynski equal- 
ity was put forward by Adib ( 2005 ) for non-Hamiltonian 



iso-energetic dynamics. It was recently generalized to en- 
ergy controlled systems by Katsuda and Ohzeki (2011 1. 



F. Weak coupling case 



In the previous sections [rvTB] [iV^Cl [RTdI we studied 
a quantum mechanical system at canonical equilibrium 
at time t = 0. During the subsequent action of the pro- 
tocol it is assumed to be completely isolated from its sur- 
rounding apart from the influence of the external work 
source and hence to undergo a unitary time evolution. 
The quality of this approximation depends on the relative 
strength of the interaction between the system and its en- 
vironment, compared to typical energies of the isolated 
system as well as on the duration of the protocol. In gen- 
eral, though, a treatment that takes into account possible 
environmental interactions is necessary. As will be shown 
below, the interaction with a thermal bath does not lead 



to a modification of the Jarzynski equality, Eq. (54), nor 
of the quantum work fluctuation relation, Eq. (53), both 



in the cases of weak and strong coupling ( Campisi et al. 



2009a Talkner et al. 2009b); a main finding which holds 
true as well for classical systems (Jarzynski 2004[ ). In 



this section we address the weak coupling case, while the 



Strictly speaking, in order to obtain well defined expressions the 
5-function has to be understood as a sharply peaked function 
with infinite support. 




FIG. 4 (Color online) Driven open system. A driven system 
[represented by its Hamiltonian 'Hs(At)] is coupled to a bath 
[represented by Hb] via the interaction Hamiltonian T-Lsb- 
The compound system is in vanishingly weak contact with a 
super-bath, that provides the initial canonical state at inverse 
temperature /3, Eq. (|60|. 



more intricate case of strong coupling is discussed in the 
next section. 

We consider a driven system described by the time de- 
pendent Hamiltonian T-lsiXt), in contact with a thermal 
bath with time independent Hamiltonian Hb, see Fig. [4] 
The Hamiltonian of the compound system is 



H{Xt)^HsiXt) + HB + HsB, 



(59) 



where the energy contribution stemming from Hsb is as- 
sumed to be much smaller than the energies of the system 
and bath resulting from 'Hs(Af) and "Hb- The parameter 
A that is manipulated according to a protocol solely en- 
ters in the system Hamiltonian, HsiXt). The compound 
system is assumed to be initially {t = 0) in the canonical 
state 



e(Ao) = e-^^<^^o)/Y(Xo), 



(60) 



where ^^(At) = Tre^'^^'-'*''^ is the corresponding partition 
function. This initial state may be provided by contact 
with a super-bath at inverse temperature /3, see Fig. |4] 
It is then assumed that either the contact is removed 
for i > or that the super-bath is so weakly coupled to 
the compound system that it bears no influence on its 
dynamics over the time span to r. 

Because the system and the environmental Hamilto- 
nians commute with each other, their energies can be 
simultaneously measured. We denote the eigenvalues of 
Hsi^t) as and those of as E^ . In analogy with 
the isolated case we assume that at time t = a joint 
measurement of HsiXo) and "Hs is performed, with out- 
comes i?^", EJ^ . A second joint measurement of 'Hs(At) 
and Hb a,t t — t yields the outcomes E^ , E^ . 

In analogy to the energy change of an isolated sys- 
tem, the differences of the eigenvalues of system and 
bath Hamiltonians yield the energy changes of system 
and bath, Ai? and AE^ , respectively, in a single realiza- 
tion of the protocol, i.e.. 



AE 
AE^ 



E^- 



(61) 
(62) 



In the weak coupling limit, the change of the energy 
content of the total system is given by the sum of the 
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energy changes of the system and bath energies apart 
from a neghgibly small contribution due to the interac- 
tion Hamiltonian Hsb- The work w performed on the 
system coincides with the change of the total energy be- 
cause the force is assumed to act only directly on the sys- 
tem. For the same reason, the change of the bath energy 
is solely due to an energy exchange with the system and 
hence, can be interpreted as negative heat, AE^ = —Q. 
Accordingly we have^^ 



the deviations being of second order in the system-bath 
interaction dTalkner et aLl|2009b|). 

Using Eqs. ( 66 1 in combination with microreversibility, 
Eq. ([io]), leads, in analogy with Eq. ([52]), 



to 



AE = w + Q. 



(63) 



Following the analogy with the isolated case, we 
consider the joint probability distribution function 
p[AE, Q; X] that the system energy changes by AE and 
the heat Q is exchanged, under the protocol A: 



ZsiXo)G[u, v; A] = Zs{\r)G[-u + il3, -v - 1/3; A] , (67) 
where 

ZsiXt) = Trse-^«-(^') , (68) 

with Tr5 denoting the trace over the system Hilbert 
space. Upon applying an inverse Fourier transform of 
Eq. ( 67 1 one arrives at the following relation: 



p[AE,Q;X] 



E 

m,n.fi.v 



5[AE - E^ + K^']<5[Q + E^ - El 



p[Ag,Q;A] 
p[-Ae-qCx] 



_ B{AE-Q-AFs) 



where 



^ Vmn\nu 



(64) 



AFs = -r'HZs{Xr)/Zs{X^)] 



(69) 



(70) 



where Pm/j|rii/[A] is the conditional probability to obtain 
the outcome E^ , Ef^ at r, provided that the outcome denotes the system free energy difference. Eq. (jeOf gen 
i?^, E^ was obtained at time t = 0, whereas p^, is eralizes the Tasaki-Crooks fluctuation theorem, Eq. ([p. 



the probability to find the outcome E^° , E^^ in the first 
measurement. The conditional probability p„^l„y [A] can 
be expressed in terms of the projectors on the common 
eigenstates of 'Hs(At), 'Hs, and the unitar y evolution gen- 
erated by the total Hamiltonian ■H(At) ( [Talkner et al. 
2009bl). 



to the case where the system can exchange heat with a 
thermal bath. 

Performing the change of the variable AE — >■ w, Eq. 
(631, in Eq. (691 leads to the following fluctuation re- 



lation for the joint probability density function of work 
and heat, namely: 



By taking the Fourier transform of p[AE, Q; X] with 
respect to both AE and Q one obtains the characteris- 
tic function of system energy change and heat exchange, 
reading 



p[w,Q;X] ^^/3(^_AFs 



p[-w,-Q-X] 



(71) 



G[u,v;X] = / d{AE)dQe 



J{uAE+vQ) 



p[AE,Q;X], (65) 



which can be further simplified and cast, as in the iso- 
lated case, in the form of a two-time quantum correlation 
function (|Talkner et aLl|2009b|): 



Notably, the right hand side does not depend on the heat 
Q but depends on the work w only. This fact implies 
that the marginal probability density of work p[w; X] = 
J dQp[w, Q; X] obeys the Tasaki-Crooks relation: 



p[w; X] 
p[-w; A] 



(72) 



G[u,v;X] = (e*["«?.(A.)-«f.le-^["«s(A„)-.«B]^^ (gg) 

where the average is over the state g{Xo), that is the 
diagonal part of g{Xo), Eq. (60), with respect to 
{'Hs(Ao), 'Hb}- Notably, in the Hmit of weak coupling 
this state g{Xo) approximately factorizes into the prod- 
uct of the equilibrium states of system and bath with 



Subsequently the Jarzynski equality. 



(531, 



is also satisfied. Thus the fluctuation relation, Eq 
and the Jarzynski equality, Eq. (54), keep holding, unal 



tered, also in the case of weak coupling. This result was 
originally found upon assuming a Markovian quantum 
dynamics for the reduced system dynamics S}^ With 



the above derivation we followed Talkner et al. (2009b) 



By use of the probability distribution in Eq. | |64| , the averaged 
quantity (AE) x = J d(AE)dQ p[AE, Q; \]AE = 1v QrHsiK)- 
Tr t?(-^o)Ws(^o) cannot, in general, be interpreted as a change 
in thermodynamic internal energy. The reason is that the final 
state, Qt, reached at the end of the protocol is typically not a 
state of thermodynamic equilibrium; hence its themrodynamic 
internal energy is not defined. 



in which one does not rely on a Markovian quantum evo- 
lution and consequently the results hold true as well for 
a general non-Markovian reduced quantum dynamics of 
the system dynamics S. 



See dCrooksl |2008| |de Roeck and Maes| |2004| [Esposito and] 
Mukamel[ |2006| |Mukamel[ [2003^ 
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G. Strang coupling case 

In the case of strong coupling, the system-bath interac- 
tion energy is non-neghgible, and therefore it is no longer 
possible to identify the heat as the energy change of the 
bath. How to define heat in a strongly coupled driven 
system and whether it is possible to define it at all cur- 
rently remain open problems. This, however does not 
hinder the possibility to prove that the work fluctuation 
relation, Eq. ( 72 ) , remains valid also in the case of strong 
coupling. For this purpose it suffices to properly iden- 
tify the work w done on, and the free energy Fs of an 
open system, without entering the issue of what heat 
means in a strong-coupling situation. As for the classical 
case, see Sec. III. A the system Hamiltonian Hs{^t) is 



the only time dependent part of the total Hamiltonian 
'H(At). Therefore the work done on the open quantum 
system, coincides with the work done on the total sys- 
tem, as in the weak coupling case treated in the previous 
section. Sec. |IV.F| Consequently, the work done on an 
open quantum system in a single realization is 



w 



E. 



(73) 



where are the eigenvalues of the total Hamiltonian 
•H(At). 

Regarding the proper identification of the free energy 
of an open quantum system, the situation is more in- 
volved because the bare partition function Z^{\t) = 
Trge"'^^-^^'*'*) cannot take into account the full effect of 
the environment in any case other than the limiting situ- 
ation of weak coupling. For strong coupling the equilib- 
rium statistical mechanical description has to be based 
on a partition function of the open quantum system that 
is given as the ratio of the partition functions of the total 
system and the isolated environment,^^ i.e.: 



(74) 



where Zb = Tr^e-'^^^ and Y{Xt) = Tre"'^^'^') with 
Trs, Tr denoting the trace over the bath Hilbert space 
and the total Hilbert space, respectively. It must be 
stressed that in general, the partition function Zs{Xt) of 
an open quantum system differs from its partition func- 
tion in absence of a bath: 



ZsiXt) + Trse- 



(75) 



The equality is restored, though, in the limit of of weak 
coupling. 

The free energy of an open quantum system follows 
according to the standard rule of equilibrium statistical 



See ijCaldeira an d Leg gett|| l983||Campisi et al. ] |2009a|b| |2010| 
Dittrich et aL[ 1998; Fcynm an and Vemon[ ^ 1963 Ford et al. 



19851 IGrabertef aL,,1988||1984nHanggi and Ingold, .2005, .2006, 



Horhammer and Buttner[|2008| l ngold[|2002||Nieuwenhui2en and| 
AUahverdyanl 120021. 



mechanics as 

Fs{Xt)^F{Xt)-FB = -lln^^. (76) 

P Zb 

In this way the influences of the bath on the thermody- 
namic properties of the system are properly taken into 
account. Besides, Eq. (FtgI) complies with all the Grand 



Laws of thermodynamics ( jCampisi et al. [2009a). 

For a total system initially prepared in the Gibbs state, 
Eq. (60 1, the Tasaki-Crooks fluctuation theorem, Eq. 



(53), applies yielding 



p[w;X] Y{Xr) 

— 6 

p[-w;X] Y{Xq) 



(77) 



Since Zb does not depend on time, the salient relation 



Y{Xr)/Y{Xo) = Zs{Xr)/Zs{Xo) 



holds, leading to: 
p[w; X] 



p[-w]X] -Zs(Ao) 



(78) 



(79) 



where /S.Fs = Fs{Xr) — Fs{Xq) is the proper free energy 
difference of an open quantum system. Since w coincides 
with the work performed on the open system the Tasaki- 
Crooks theorem, Eq. (72), is recovered in the case of 



strong coupling ( Campisi et al. . 2009a ) . 



V. QUANTUM EXCHANGE FLUCTUATION RELATIONS 

The transport of energy and matter between two reser- 
voirs that stay at different temperatures and chemical 
potentials represents an important experimental set-up, 
also Sec. VI below, as well as a central prob- 



see 



lem of non-equilibrium thermodynamics (de Groot and 
Mazur 1984). Here, the two-measurement scheme de- 



scribed above in conjunction with the principle of mi- 
croreversibility leads to fluctuation relations similar to 
the Tasaki-Crooks relation, Eq. (53), for the proba- 
bilities of energy and matter exchanges. The resulting 
fluctuation relations have been referred to as "exchange 
fluctuation theorems" (Jarzynski and Wojcik 2004), to 
distinguish them from the "work fluctuation theorems" . 

The first quantum exchange fluctuation theorem was 
put forward by 'Jarzynski and Wojcik ( 2004[ ). It applies 
to two systems initially at different temperatures that 
are allowed to interact over the lapse of time (0, r), via a 
possibly time-dependent interaction. This situation was 
later general i zed b y Saito and Utsumi ( [2008) and An- 
drieux et al. (2009), to allow for the exchange of energy 



and particles between several interacting systems initially 
at different temperatures and chemical potentials; see the 
schematic illustration in Fig. [5] The total Hamiltonian 
'H(Vt) consisting of s subsystems is: 



(80) 
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FIG. 5 (Color online) Exchange fluctuation relation set-up. 
Several reservoirs (large semicircles) interact via the coupling 
Vt (symbolized by the small circle), which is switched on at 
time t = and switched oflF at t = r. During the on-period 
(0, t) the reservoirs exchange energy and matter with each 
other. The resulting net energy change of the i-th reservoir 
is AEi and its particle content changes by ANi. Initially 
the reservoirs have prescribed temperatures, Ti — {ksPi)'^^, 
and chemical potentials, fit, of the particle species that are 
exchanged, respectively. 



where Hi is the Hamiltonian of the i-th systeni, and 
Vf describes the interaction between the subsystems, 
which sets in at time t = and ends at time t = t. 
Consequently Vf = for f ^ £^nd in particular 

Vq — Vr — 0. As before, it is very important to distin- 
guish between the values Vt at a specific time and the 
whole protocol V. 

Initially, the subsystems are supposed to be isolated 
from each other and to stay in a factorized grand canon- 
ical state 



Qo = 



(81) 



with fj,i, Pi, Si — Tji^g-PiCH-i-t^i^i) the chemical potential, 
inverse temperature, and grand potential, respectively, of 
subsystem i. Here, A/i and Tr^ denote the particle num- 
ber operator and the trace of the i"^ subsystem, respec- 
tively. 

We also assume that in absence of interaction the 
particle numbers in each subsystem are conserved, i.e., 
[Hi,Afi] — 0. Since operators acting on Hilbert spaces 
of different subsystems commute, we find [Hi^Afj] = 0, 
[A/i,A/'j] = 0, and ['Hi,'Hj] = for any Accordingly, 
one may measure all the Hi^s and all the A/i's simultane- 
ously. Adopting the two-measurement scheme discussed 
above in the context the work fluctuation relation, we 
make a first measurement of all the Hi 's and all the Mi 's 
at i = 0. Accordingly, the wave function collapses onto 
a common eigenstate \tpn) of all these observable with 
eigenvalues i?^, iV,*. Subsequently, this wave function 
evolves according to the evolution ?7t,o[V] generated by 
the total Hamiltonian, until time t when a second mea- 
surement of all "Hi's and A/i's is performed leading to 
a wave function collapse onto an eigenstate lipm), with 



eigenvalues E^^, N^. As in the case studied in Sec. |IV.F 
the joint probability density of energy and particle ex- 
changes p[AE, AN; V] completely describes the effect of 
the interaction protocol V 

p[AE, AN; V] = ^ n ^(^^^ ' ^™ + 

m,n i 

x6{AN, iV;)p„|„[VK , (82) 

where Pmlni"!^] is the transition probability from state 

to 1^™) and pI = H,e-ft[^"-''-^;VS. is the ini- 
tial distribution of energies and particles. Here the sym- 
bols AE and AN, are short notations for the individ- 
ual energy and particle number changes of all subsys- 
tems ASi, AE'a, . . . , AEs and AiVi, AiVa, . . . , AA^^, re- 
spectively. 

Assuming that the total Hamiltonian commutes with 
the time reversal operator at any instant of time, and 
using the time reversal property of the transition proba- 
bilities, Eq. ([56|), one obtains that 



p[AE, AN;V] 
p[-AE,-AN;V] 



,l3i[AEi-p.iANi] 



(83) 



This equation was derived by Andrieux et al. (2009), and 



expresses the exchange fluctuation relation for the case 
of transport of energy and matter. 

In the case of a single isolated system (s 1), it re- 
duces to the Tasaki-Crooks work fluctuation theorem, 
Eq. (53), upon rewriting AEi = w and assuming that 
the total number of particles is conserved also when the 
interaction is switched on. i.e., [H{Vt),JV] — 0, to obtain 
AA^ = 0. The free energy difference does not appear in 
Eq. ( 83 ) because we have assumed the protocol V to be 



cychc. 

In the case of two weakly interacting systems {s — 2,V 
small), that do not ex change particles, it reduces to the 
fluctuation theorem of Jarzynski and Wojcik (2004) for 
heat exchange: 



p[Q;Vl 



(84) 



where Q — AEi = —AE2, with the second equality fol- 
lowing from the assumed weak interaction. 

In case of two weakly interacting systems {s — 2, V 
small) that do exchange particles, and are initially at 
same temperature, yielding Q ~ 0, the fluctuation rela- 
tion takes on the form: 



p[q;Vl 

p[-q;V] 



— g/3(Aii-P2)g 



(85) 



where q = ANi = -AN2. 

One basic assumption leading to the exchange fluc- 
tuation relation, Eq. (831, is that the initial state is 



a factorized state, in which the various subsystems are 
uncorrelated from each other. In most experimental sit- 
uations, however, unavoidable interactions between the 
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systems, would lead to some correlations, and a conse- 
quent deviation from the assumed factorized state, Eq. 



(81 1. The resulting deviation from the exchange fluctua- 



tion relation, Eq. (83 1, is expected to vanish for observa- 
tion times T larger than some characteristic time scale Tc, 
determined by the specific physical properties of the ex- 



perimcntal set-up (Andrieux et al. 
[2010a, Esposito e^ a/.[ 120091 1 



2009 Campisi et al. 



p[AE, AN; V] 
p[-AE,-AN;V] 



,pi{AEi-fiiANi) 



(86) 



For those large times t ^ Tc a non equilibrium steady 
state sets in under the condition that the reservoirs are 



chosen macroscopic. For this reason Eq. (86) is referred 
to as a "steady state fluctuation relation" . This is in con- 
trast to the other fluctuation relations discussed above, 
which instead are valid for any observation time r, and 
are accordingly referred to as transient fluctuation rela- 
tions. Saito and Dhar (2007) provided an explicit demon- 



stration of Eq. ( 86 ) for the quantum heat transfer across 
a harmonic chain connecting two thermal reservoirs at 
different temperatures. Ren et al. (20101 reported on 



the breakdown of Eq. ( 86 ) induced by a nonvanishing 



Berry-phase heat pumping. The latter occurs when the 
temperatures of the two baths are adiabatically modu- 
lated in time. 



VI. EXPERIMENTS 

A. Work fluctuation relations 

Regarding the experimental validation of the work fluc- 
tuation relation, a fundamental difference exists between 
the classical and the quantum regime. In classical ex- 
periments work is accessible by measuring the trajectory 
Xt of the possibly open system, and integrating the in- 
stantaneous power according to = / dtdHg /dt, Eq. 
(31). In clear contrast, in quantum mechanics the work 



is obtained as the difference of two measurements of the 
energy, and an "integrated power" expression does not 
exists for the work, see Sec. |III.C| 

Following closely the prescriptions of the theory one 
should perform the following steps in order to experi- 
mentally verify the work fluctuation relation, Eq. (53): 



(i) Prepare a quantum system in the canonical state, Eq. 
([2]), at time t — 0. (ii) Measure the energy at t — 0. (hi) 
Drive the system by means of some forcing protocol At 
for times t between and r, and make sure that during 
this time the system is well insulated from its environ- 
ment, (iv) Measure the energy again at r and record 
the work w, according to Eq. (39). (v) Repeat this pro- 
cedure many times and construct the histogram of the 
statistics of work as an estimate of the work pdf p[w; A]. 
In order to determine the backward probability the same 
type of experiments has to be repeated with the inverted 
protocol, starting from an equilibrium state at inverse 



temperature /3 and at those parameter values that are 
reached at the end of the forward protocol. 



1. Proposal for an experiment employing trapped cold ions 



Huber et al. { 2008 1 suggest an experiment that follows 
exactly the procedure described above. They propose to 
implement a quantum harmonic oscillator by optically 
trapping an ion in the quadratic potential generated by 
a laser trap, using the set-up devel oped by|Schulz et al] 
(|2008D. In principle, the set-up of [Schulz et al.\ (|2008[) 



allows, on one hand, to drive the system by changing in 
time the stiffness of the trap, and, on the other hand, to 
probe whether the ion is in a certain Fock state |n) ; i.e., in 
an energy eigenstate of the harmonic oscillator. The mea- 
surement apparatus may be understood as a single Fock 
state "filter" , whose outcome is "yes" or "no" , depending 
on whether the ion is or is not in the probed state. Thus 
the experimentalist probes each possible outcome (n, m), 
where (n, to) denotes the Fock states at time t = and 
t = T, respectively. Then, the relative frequency of the 
outcome (n, m) occurrence is recored by repeating the 
driving protocol many times always preparing the sys- 
tem in the same canonical initial state. In this way the 
joint probabilities Pm|n[A]p^ are measured. 

The work histogram is then constructed as an estimate 
of the work pdf, Eq. ( 48 ) , thus providing experimen- 

Like- 



(531 



tal access to the fluctuation relation, Eq. 
wise the relative frequency of the outcomes having n 
as the initial state gives the experimental initial popu- 
lation p° . Thus, with this experiment one can actually 
check the symmetry relation of the conditional proba- 
bilities, Pm|n[A] = Pn|m[A], Eq. (56) and compare their 



experimental values with the known theoretical values 



(Deffner and Lutz 


2008a 


2009a 


)■ 



2008[ |Husimi[ |1953[ [Talkner et af. 



Another suitable quantum system to test quantum 
fluctuation relations are quantum versions of nanome- 
chanical oscillator set-ups that with present day nan- 
otechnology are at the verge of entering the quantum 
regime. In these systems work protocols can be im- 
posed by optomechanical means. 



2. Proposal for an experiment employing circuit quantum 
electrodynamics 



Currently, the experiment proposed by |Huber et al.\ 
(2008) has not yet been carried out. An analogous exper- 



iment could, in principle, be performed in a circuit quan- 
tum electrodynamics (QED) set-up as the one described 
in (IHofheinz et aL||2008l 120091). The set-up consists of 



Note the exciting recent advancements obtained with the works 
( Anctsbcr ger et aL| [20091 [Kippenberg and Vahala| |2008| |LaHaye| 
I et a;.||2004||O'Connell ef aL[|2010[ |. 
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a Cooper pair box qubit (a two states quantum system) 
that can be coupled to and de-coupled from a supercon- 
ducting ID transmission line, where the latter mimics a 
quantum harmonic oscillator. With this architecture it 
is possible to implement various functions with very high 
degree of accuracy. Among them the following tasks are 
of special interest in the present context: (i) Creation of 
pure Fock states i.e., the energy eigenstates of the 
quantum harmonic oscillator in the resonator, (ii) Mea- 
surement of photon statistics Pm, i-e., measurements of 
the population of each quantum state \m) of the oscilla- 
tor, (iii) Driving the resonator by means of an external 
field. 



Hofheinz et al. (2008) report, for example, on the cre- 
ation of the ground Fock state |0), followed by a driv- 
ing protocol A (a properly engineered microwave pulse 
applied to the resonator) that "displaces" the oscillator 
and creates a coherent state |a), whose photon statis- 

to 



tics Pm|o[-'^] measured with good accuracy up 
n-max ^ 10. In more recent experiments (Hofheinz et al. 



2009 ) the accuracy was improved and n^ax has raised to 
~ 15. The quantity Pm|o[-^] is actually the conditional 
probability to find the state |m) at time t — t, given 
that the system was in the state |0) at time t — 0. Thus, 
by preparing the oscillator in the Fock state \n) instead 
of the ground state |0), and repeating the same driving 
and readout as before, the matrix Pm|„[A] can be de- 
termined experimentally. Accordingly one can test the 
validity of the symmetry relation Pm|n[A] = PTt|m[A], Eq. 
(56), which in turn implies the work fluctuation relation, 
see Sec. |IV.E[ At variance with the proposal of |Huber] 
et al. (2008), in this case the initial state would not be 



randomly sampled from a canonical state, but would be 
rather created deterministically by the experimenter. 



The theoretical values of transition probabilities for 
this case corresponding to a displacement of the oscillator 



were first reported by ( Husimi 1953), see also (Campisi 



2008). Talkner et al. (2008a) provided an analytical ex- 
pression for the characteristic function of work and inves- 
tigated in detail the work probability distribution func- 
tion and its dependence on the initial state, such as for 
example canonical, microcanonical, and coherent state. 
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FIG. 6 (Color online) Scheme of a bidirectional counting 
statistics experiment. Two leads (large semicircles) with dif- 
ferent electronic chemical potentials (/ii ^ (12) and same tem- 
perature (/3i =132) are connected through a double quantum 
dot (small circles), whose quantum state is continuously mon- 
itored. The state (1,0), i.e., "one electron in the left dot, no 
electrons in the right dot" is depicted. The transition from 
this state to the state (0,1) signals the exchange of one elec- 
tron from subsystem 1 to subsystem 2. "Hi. 2 and A/1,2 denote 
the Hamiltonian and electron number operators of the sub- 
systems, respectively. 



B. Exchange fluctuation relations 

Like the quantum work fluctuation relations, the quan- 
tum exchange fluctuation relations are understood in 
terms of two point quantum measurements. In an ex- 
perimental test, the net amount of energy and/or num- 
ber of particles (depending which of the three exchange 
fluctuation relations, 



Eqs. (|83) [84) [85]), is studied) has 
to be measured in each subsystems twice, at the be- 
ginning and at the end of the protocol. However typi- 
cally these are macroscopic reservoirs, whose energy and 
particle number measurement is practically impossible. 
Thus, seemingly, the verification of the exchange fluctu- 
ation relations would be even more problematic than the 
validation of the quantum work fluctuation relations. In- 
deed, while experimental tests of the work fluctuation re- 
lations have not been reported yet, experiments concern- 
ing quantum exchange fluctuation relations have been 
already performed. In the following we shall discuss two 
Utsumi et al.\ ( [2010 ) and the other by 



of them, one by 



Nakamura et al. (2010). In doing so we demonstrate how 
the obstacle of energy/particle content measurement of 
macroscopic reservoirs was circumvented. 



1. An electron counting statistics experiment 



So far we have addressed possible experimental tests of 
the Tasaki-Crooks work fluctuation theorem, Eq. (53), 



for isolated systems. The case of open systems, interact- 
ing with a thermal bath, poses extra difficulties related 
to the fact that in order to measure the work in this case 
one should make two measurements of the energy of the 
total macroscopic system, made up of the system of inter- 
est and its environment. This presents an extra obstacle 
that at the moment seems difficult to surmount except 
for a situation at (i) weak coupling and (ii) Q ^ 0, then 
yielding, together with Eq. (63) w ^ AE. 



Utsumi et al. (2010) have recently performed an exper- 
imental verification of the particle exchange fluctuation 
relation, Eq. ( 85 ) , using bi-directional electron counting 
statistics (Fujisawa et al. 2006). The experimental set- 



up consists of two electron reservoirs (leads) at the same 
temperature. The two leads are connected via a double 
quantum dot, see Fig. [6] When an electric potential dif- 



See (Andrieux et a/. "2009 !|Campisi et aI]|2010a|[Esposit o et al 



20091, and also |JarzyiiskiJ 2000 1 regardini 



the classical case. 
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ference V = fii — fi2 applied to the leads, a net flow 
of electrons starts transporting charges from one lead to 
the other, via lead-dot and dot-dot quantum tunnelings. 
The measurement apparatus consists of a secondary cir- 
cuit in which a current flows due to an applied voltage. 
Thanks to a properly engineered coupling between sec- 
ondary circuit and the double quantum dot, the current 
in the circuit depends on the quantum state of the dou- 
ble dot. The latter has four relevant states, which we 
shall denote as |00), |01), |10), |11) corresponding respec- 
tively to: no electrons in the left dot and no electrons in 
the right dot, no electrons in the left dot and one elec- 
tron in the right dot, etc.,. Each of these states leads to 
a different value of the current in the secondary circuit. 
In the experiment an electric potential difference is ap- 
plied to the two leads for a time t. During this time the 
state of the double quantum dot is monitored by regis- 
tering the current in the secondary circuit. This current 
was found to switch between the four values correspond- 
ing to the four quantum states mentioned above. The 
outcome of the experiment is a sequence Ik of current 
values, with 1^ taking only four possible values. In other 
terms, the outcome of the experiment consists of a se- 
quence {{l,r)}k (with l,r = 0,1) of joint eigenvalues of 
two commuting observables £,TZ specifying the occupa- 
tion of the left (l) and right (r) dots by single electrons 
at the time of the fc"^ measurement. The presence of 
an exchange of entries within one time step of the form 
(1, 0)„, (0, l)ii+i signals the transfer of one electron from 
left to right, and vice versa (0, 1)„, (l,0)„+i the transfer 
from right to left. Thus, given a sequence {{l,r)}k, the 
total number q[{{l, r)}k] of electrons transferred from left 
to right is obtained by subtracting the total number of 
right-to-left transfers from the total number of left-to- 
right transfers. It was found that, for observation times 
larger than a characteristic time Tc, the fluctuation rela- 



tion, p[q] ~ p[—q\e"^'^, Eq. (85), was satisfied with the 



actual temperature of the leads replaced by an effective 
temperature, see Fig [T] The renormalization of temper- 
ature was explained as an effect due to an exchange of 
electrons occurring between the dots and the secondary 



circuit (Utsumi et at. 2010). 



The question however remains of how to connect this 
experiment in which the flux of electrons through an in- 
terface is monitored and the theory, leading to Eq. ( 85 ) , 



which instead prescribes only two measurements of total 
particle numbers in the reservoirs. The answer was given 
in Campisi et al. (2010a), who showed that the exchange 
fluctuation relation, Eq. (83), remains valid, if in addi- 



tion to the two measurements of total energy and particle 
numbers occurring at and t, the evolution of a quan- 
tum system is interrupted by means of projective quan- 
tum measurements of any observable A that commutes 
with the quantum time reversal operator Q. In other 
words, while the forward and backward probabilities are 
affected by the occurrence of intermediate measurement 
processes, their ratio remains unaltered. 



need to measure the initial and final content of particles 
in the reservoirs because the number of exchanged par- 
ticles is inferred from the sequence of intermediate mea- 
surements outcomes {(Z,r)}fc. Thus, thanks to the fact 
that quantum measurements do not alter the fluctuation 
relation, one may overcome the problem of measuring 
the energy and number of particles of the macroscopic 
reservoirs, by monitoring instead the flux through a mi- 
croscopic junction. 



2. Nonlinear response relations in a quantum coherent 
conductor 

As discussed in the introduction, the original motiva- 
tion for the study of fluctuation relations was to overcome 
the limitations of linear response theory and to obtain 
relations connecting higher order response functions to 
fluctuation properties of the unperturbed system. As an 
indirect and partial confirmation of the fluctuation re- 
lations higher order static fluctuation-response relations 
can be tested experimentally. 

Such a validation was recently accomplished in co- 
herent quantum transport experiments by |Nakamura| 



et al. (2010 2011), where the average current / and 



the zero-frequency current noise power S generated in 
an Aharonov-Bohm ring were investigated as a function 
of an applied dc voltage V , and magnetic fleld B. In the 
nonlinear response regime, the current and noise power 
may be expressed as power series of the applied voltage: 



J(l/,B) = Gi(B)l/- 



G2(B)^2 , G3(B) 



3! 



(87) 



S{V, B) = 5o(B) + 5i(B)y + ^^V^ 



where the coefficients depend on the applied magnetic 
field B. The steady state fluctuation theorem, Eq. (86), 



then predicts the following fluctuation relations (Saito 
and Utsumil 120081) 



So = ^kBTGi, Sf^2kBTGi, Sf ^ GkeTG^ (89) 

where 5f = S'i(B) + S',(-B), Sf = S',(B) - S',(-B), and 
analogous definitions for Gf and Gf. T he first equation 



in (89) is the Johnson-Nyquist relation (Johnson 1928 
Nyquist 1928). In the experiment by Nakamura et al. 



(2010) good quantitative agreement with the first and 



the third expressions in (89) was established, whereas. 



for the time being, only qualitative agreement was found 
with the second relation. 

The higher order static fiuctuation dissipation rela- 
tions ( 89 ) were obtained from a steady state fluctuation 



theorem for particle exchange under the simplifying as- 



sumption that no heat exchange occurs ( Nakamura et al. 
2010|. Then the probability of transferring q particles 



In the experiment of Utsumi et al. (2010) one does not 



is related to the probability of the reverse transfer by 
p{q) = p{—q)e^'^ where A = (iV = /3(/ii — ^2) is the 
so-called affinity. If both sides are multiplied by q and 
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integrated over q a comparison of equal powers of applied 



voltage y yields Eq. ( 89 ) ( Nakamura ei a/. 2011). An al- 
ternative approach, that also allows to include the effect 
of heat conduction, is offered by the fluctuation theorems 
for currents in open quantum systems. This objective has 
be en put forward by Saito and Utsumi (2008) and also 
by Andrieux et al. f 2009) , based on a generating function 
approach in the spirit of Eq. ( 55 ) . 



VII. OUTLOOK 

In closing this colloquium we stress that the known 
fluctuation relations are based on two facts: (a) microre- 
versibility for non-autonomous Hamiltonian systems Eq. 
(40), and (b) the special nature of the initial equilib- 
rium states which is expressible in either micro-canonical, 
canonical or grand-canonical form, or products thereof. 
The final state reached at the end of a protocol though 
is in no way restricted. It evolves from the initial state 
according to the governing dynamical laws under a pre- 
scribed protocol. In general this final state may markedly 
differ from any kind of equilibrium state. 

For quantum mechanical systems it also is of utmost 
importance to correctly identify the work performed on 
a system as the difference between the energy of the sys- 
tem at the end and the beginning of the protocol. In case 
of open systems the difference of the energies of the total 
system at the end and beginning of the protocol coincides 
with the work done on the open system as long as the 
forces exclusively act on this open system. With the free 
energy of an open system determined as the difference 
of free energies of the total system and that of the iso- 
lated environment the quantum and classical Jarzynski 
equality and the Tasaki-Crooks theorem continue to hold 
true even for systems strongly interacting with their envi- 
ronment. Deviations from the fluctuation relations how- 
ever must be expected if protocol forces not only act on 
the system alone but as well directly on the environmen- 
tal degrees of freedom; for example, if a time-dependent 
system-bath interaction protocol is applied. 

The most general and compact formulation of quantum 
work fluctuation relations also containing the Onsager- 
Casimir reciprocity relations and nonlinear response to 
all orders, is the Andrieux-Gaspard relation, Eq. (55) 



which represents the proper quantum version of the clas- 



sical Bochkov-Kuzovlev formula ( Bochkov and Kuzovlev 
1977), Eq. (14). These relations provide a complete 



theoretical understanding of those nonequilibrium situ- 
ations that emerge from arbitrary time-dependent per- 
turbations of equilibrium initial states. 

Less understood are exchange fluctuation relations 
with their important applications to counting statistics 
(Esposito et al. 2009). The theory there so far is re- 



stricted to situations where the initial state factorizes into 
grand-canonical states of reservoirs at different temper- 
atures or chemical potentials. The interaction between 
these reservoirs is turned on and it is assumed that it will 



lead to a steady state within the duration of the protocol. 
Experimentally, it is in general difficult to exactly follow 
this prescription and therefore a comparison of theory 
and experiment is only meaningful for the steady state. 
Alternative derivations of exchange relations for more re- 
alistic, non-factorizing initial states would certainly be of 
much interest. In this context, the issue of deriving quan- 
tum fluctuation relations for open systems that initially 
are in nonequilibrium steady quantum transport states 
constitutes a most interesting challenge. Likewise, from 
the theoretical point of view little is known thus far about 
quantum effects for transport in presence of time depen- 
dent reservoirs, for example using a varying temperature 
and/or chemical potentials (Ren et al. 2010). 



The experimental applications and validation schemes 
involving nonlinear quantum fluctuation relations still 
are in a state of infancy, as detailed with Sec. |VI[ so 
that there is plenty of room for advancements. The ma- 
jor obstacle for the experimental verification of the work 
fluctuation relation is posed by the necessity of perform- 
ing quantum pr ojective measuremen ts of energy. Besides 
the proposal of Huber et al. ( |2008 ) employing trapped 
ions, we suggested here the scheme of a possible exper- 
iment employing circuit-QED architectures. In regard 
to exchange fluctuation relations instead, the main prob- 
lem is related to the difhculty of measuring microscopic 
changes of macroscopic quantities pertaining to heat and 
matter reservoirs. Continuous measurements of fluxes 
seemingly provide a practical and efhcient loophole for 
this dilemma (Campisi et al. 2010a). 



The idea t hat useful work may be obtained by using 
information ( Maruyama et al. . 2009 ) has established a 



connection between the topical fields of quantum infor- 
mation th eory dVedral, ,2002) and quantum fluctuation 



relations. Piechocinska (2000) and Kawai et al. (2007) 



used fluctuation relations and information theoretic mea- 
sures to derive Landauer's principle. A generalization of 
the Jarzynski equality to the case of feedback controlled 



systems was provided in the classical case by'Sagawa and 



Ucda (20 10l, and in the quantum case by .Mo rikuni and 
[Tasalkil f 20111 . Recently [Peffner and Lut'zn2010[ ) gave 



bounds on the entropy production in terms of quantum 
information concepts. In similar spirit, |Hide and Vedral] 
(2010) presented a method by relating relative quantum 



entropy to the quantum Jarzynski fluctuation identity in 
order to quantify multi-partite entanglement within dif- 
ferent thermal quantum states. A practical application 
of the Jarzynski equality in quantum computation was 
showed by lOhzekil (|2010| . 



In conclusion, the authors are confident in their belief 
that this topic of quantum fluctuation relations will ex- 
hibit an ever growing activity within nanosciences and 
further may invigorate readers to pursue still own re- 
search and experiments as this theme certainly offers 
many more surprises and unforeseen applications. 
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Appendix A: Derivation of the Bochkov-Kuzvlev relation 



We report below the steps leading to Eq. ( 14 ) 



exp 



dzo- 







Z{to) 
dzT-po(zr)cxp 

dZrPo(Zr)cxp 



exp 



ds UsB{ips,o[^o; A]) 



T 

dr Ur-r£BB{iprfi[z'^; £qX]) 



(Al) 



where the first equality provides an explicit expression 
for the l.h.s. of Eq. (14 1. In going from the second to 



the third line we employed the expression of work in Eq. 
(131, the microreversibility principle (10) and made the 
change of variable Zq z.^- The Jacobian of this trans- 
formation is unity, because the time evolution in classi- 
cal mechanics is a canonical transformation. A further 
change of variables z^- — > z'^ = ez^-, whose Jacobian is 
unity as well, and the change s ^ r — t — s, yields the 
expression in the last line, that coincides with the right 
hand side of Eq. (14). In the last line we used the prop- 



erty po(z) = pais^Y, inherited by po = e /Zq from 
the assumed time reversal invariance of the Hamiltonian, 
H{z) = H{ez). 



Appendix B: Quantum microreversibility 

In order to prove the quantum principle of microre- 
versibility, we first discretize time and express the time 
evolution opera tor Ut.o[X] as a time ordered product 
( |Schleich[ [200T| ) 



N—>-oo 

(Bl) 

where e = t/N denotes the time step. Using Eq. (|9|, we 
obtain 



C/r-t,o[A] = lim e 



(B2) 



Thus: 

e^Ur-tM& = (B3) 



where we inserted 80^ = 1, N — 1 times. Assuming 
that -ff (At) commutes at all times with the time reversal 
operator 6, Eq. (|4T|), we find 



(B4) 



for any complex number u. Using this equation and the 
fact that £ is real- valued, e* = e, we obtain Eq. (40): 



lim 



t 



C/;jA] = C/t,.[A]. 



(B5) 



Appendix C: Tasaki-Crooks relation for the characteristic 
function 



From Eq. ( 50 ) we have 



Z{Xq)G[u;X] = TrC/^o[A]e*"^(^-)t/^o[A]e-"'^(^°)e-'^«(^«). 

(CI) 



For t — 0, the microreversibility principle, Eq. (40), be- 



comes Uo^r[X] = C/]_o[A] = Q^Ur,o[X]&. Therefore, 
Z{Xo)G[u;X] = 

Tr Ur [A] ee™^(^- ^ 6^ C/^ q [A] Oe-^^^^") e-''^^^^^ 9^ 6 , 

(C2) 

where we inserted Q^Q = 1 under the trace. Using Eq. 
(B4), we obtain 



Z{Xo)G[u:X]= (C3) 
Tre^t/,^o[A]e-™'^(^-)?7^_o[A]e™'^^^°^e-'3^(^«)e. 

The anti-linearity of O implies, for any trace class oper- 
ator A: TrOMe ^TtA^ . Using this we can write 

Z{Xo)G[u;X] = Tre-'^^(^«)e-™^(^°)[/,o[A]e™^(^-)t/^o[A] 

(C4) 

Using the cyclic property of the trace one then obtains 
the important result 

Z{Xo)G[u;X] 

= Tr[/|o[A]e*(-"+^««(^«)C/,,o[A]e-'(-"+^^)^(^°)e-^«(^'') 
^Z{Xr)G[-u + i(5;X] . (C5) 
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Addendum and Erratum: Colloquium: Quantum Fluctuation Relations: 
Foundations and Applications [Rev. Mod. Phys. 83, 779 (2011)] 

Michele Campisi, Peter Hanggi, and Peter Talkner 

Institute of Physics, University of Augsburg, Universitatsstr. 1, D-86135 Augsburg, Germany 
(Dated: July 26, 2011) 



PACS numbers: 05.30.-d, 05.40.-a 05.60.Gg 05.70.Ln, 



The first line of Eq. (51) contains some typos: it correctly reads 

G[u: A] = Trre™[«?(^-'-«(^'')le-''«(^''V-2(Ao) • (51) 



This compares with its classical analogue, i.e, the second line of Eq. (27). 
Quite surprisingly, notwithstanding the identity 



Jo 



dXt 

one finds that generally 

dH?(Xt 



L ^0 "^t 



(a) 



(b) 



As a consequence it is not allowed to replace ^^{Xr) — ^{{Xq), with AtXtdT-L^ (Xt) / dXt in Ecj. (51). Tims, there 
is no quantum analogue of the classical expression in the third line of Eq. (27). This is yet another indication that 
"work is not an observable" (Talkner et at, 2007). This observation also corrects the second line of Eq. (4) of the 
original reference (Talkner et al, 2007). 



The correct expression is obtained from the general formula 



Texp [A{t) - A{0)] = Texp 



-'0 



at 



(c) 



where A(t) is any time dependent operator [in our case A{t) = iu'Hf{Xt)]- Equation (c) can be proved by demonstrat- 
ing that the operator expressions on either side of Eq. (c) obey the same differential equation with the identity operator 
as the initial condition. This can be accomplished by using the operator identity de^^*^ /dt = Jq ds e^^^*^ A{t)e^^~^^^^*\ 

There are also a few minor misprints: (i) The symbol ds in the integral appearing in the first line of Eq. (55) 
should read dt, (ii) the correct year of Reference (Morikuni and Tasaki, 2010) is 2011 (not 2010). 

The authors are grateful to Prof. Yu. E. Kuzovlev for providing them with this insight, and for pointing out the 
error in the second line of Eq. (51). 
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